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ABSTRACT 

Compact group galaxies often appear unaffected by their unusually dense environment. Closer 
examination can, however, reveal the subtle, cumulative effects of multiple galaxy interactions. 
Hickson Compact Group (HCG) 59 is an excellent example of this situation. We present a 
photometric study of this group in the optical (HST), infrared (Spitzer) and X-ray (Chandra) 
regimes aimed at characterizing the star formation and nuclear activity in its constituent galaxies 
and intra-group medium. We associate five dwarf galaxies with the group and update the velocity 
dispersion, leading to an increase in the dynamical mass of the group of up to a factor of 10 (to 
2.8 x 10 13 Mq), and a subsequent revision of its evolutionary stage. Star formation is proceeding 
at a level consistent with the morphological types of the four main galaxies, of which two are 
star-forming and the other two quiescent. Unlike in some other compact groups, star-forming 
complexes across HCG 59 closely follow mass-radius scaling relations typical of nearby galaxies. 
In contrast, the ancient globular cluster populations in galaxies HCG 59A and B show intriguing 
irregularities, and two extragalactic H II regions are found just west of B. We age-date a faint 
stellar stream in the intra-group medium at ~ 1 Gyr to examine recent interactions. We detect a 
likely low-luminosity AGN in HCG 59A by its ~ 10 40 erg s" 1 X -ray emission; the active nucleus 
rather than star formation can account for the UV+IR SED. We discuss the implications of our 
findings in the context of galaxy evolution in dense environments. 

Subject headings: galaxies: clusters: individual (HCG 59) — galaxies: star clusters — galaxies: evolution 
— galaxies: interactions — galaxies: active — galaxies: dwarf — galaxies: fundamental parameters 
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1. Introduction 

Compact galaxy groups populate the high den- 
sity tail of the galaxy number density distribution. 
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The systems catalogued by IfficksoiJ (|l982i Hick- 
son Compact Groups, or HCGs) exhibit some fea- 
tures, such as dynamical and evolutionary states, 
elliptical fractions and X-ray properties of the 
intra-group medium (IGM) similar to galaxy clus- 
ters, their massive, more populous counterparts. 
In contrast to the well-studied cluster galaxies, 
however, the specific effects of the compact group 
environment on the evolution of its galaxies are 
not yet clear. 

HCGs are defined through criteria of isola- 
tion and surface brightnes^l that give rise to self- 
gravitating, dense groupings of a few (typically 
four) main members. Because of their masses, 
these galaxies orbit around the group barycen- 
ter rather sluggishly, with velocity dispersions 
on the order of <jcg ~ 250 kms -1 (jTago et al 



20081 : ICoxl l2000h . cf. galax y cluster dispersions 



of fJ, 



Hnst.fir 



750 kms" 1 (IBinggeli et al.1 11987 



The fc White! Il986l [Cox|l2000|). This trait makes 



HCGs valuable laboratories for galaxy evolution: 
the low velocity dispersions force some galax- 
ies into strong, prolonged interactions while oth- 
ers appear undisturbed but are apparently un- 
dergoing enhanced secular evolution. That is to 
say, this latter population is affected by gravita- 
tional interplay with their neighbors, but evolve 
more subtly, without obvious , strong interactions 
( Konstantopoulos et al. 2010h . 

Relating the various observational characteris- 
tics of compact groups to those of clusters is im- 
portant for understanding whether they constitute 
their own class, or if they are simply mini-clusters. 
Perhaps more appropriately, structures like com- 
pact groups may be considered pla usible buildin, 
blocks of clusters at high er z (e. g. iFuiita fc Got 



2004 iRudick et al.l 120061 ). Revealing past inves 



tigations of HCGs as a class have focussed on 
gas content. Their members are typically defi- 
cient in H I gas when compared to galaxies of 
similar morphological types and masses (e. g. the 
samp le of isolated galaxies in Haynes fc Giovanelli 



19841 ) . IVerdes-Montenegro et al.l (|200ll ) proposed 



an evolutionary sequence based on mapping the 
spatial distribution of Hi ac ross a large sample of 
HCGs. I Johnson et"aT1 (j2007l hereafter J07) added 



1 $iv > 3 9q, i. e. a circular area defined by three galaxy- 
mean-radii about the group is devoid of galaxies of com- 
parable brightness. A group surface brightness of fj, < 
26.0 mag defines galaxy density. 



to this investigation by quantifying the gas rich- 
ness of twelve groups with the relation of H I- 
to-dynamical mass, log(Aim)/ log(Afdyn)- This 
gave rise to the hypothesis of an alternate, two- 
pronged ev olutionary diagram for HCGs, w hich we 
explored in Konstantopoulos et al. ( 2010l) . In one 
path, the galaxies have strong interactions before 
exhausting their cold gas reservoirs for star forma- 
tion, in the other, gas is processed by star forma- 
tion within individual galaxies prior to late-stage 
dry mergers. 

Furthermore, the mid-IR colors of HCG galax- 
ies show an interesting bimodal distribution that 
distinguishes star- f orming from quiescent systems. 
I Walker et al.l (|2010l ) interpret this statistically sig- 
nificant gap as evidence for accelerated galaxy evo- 
lution in the compact group environment. Their 
similar mid-IR color distributions relate HCGs to 
the infall regions of clusters and set them apart 
from any other galaxy sample compared, interact- 
ing or quiescent. This theme was expanded by 
Tzanavaris et al.l (|2010h who found this gap appar- 



ent also in the distribution of specific star forma- 
tion rates for HCG galaxies. These observations 
together point to compact groups as local exam- 
ples of the plausible building blocks of clusters in 
the early universe. 

In addition, HCGs, which are isolated by selec- 
tion, could potentially help explain the evolution- 
ary history of some field ellipticals. For example, 
Rubin et al.l ( 1990T) o riginally proposed (see also 



Gallaghe r et al.l 120081 ) that HCG 31 will evolve 



into a single, field elliptical through a wet merger 
(one where gas is still available during the interac- 
tion). 'Fossil groups', the probable ulti mate fate of 
i solate d groupings, were examined by I Jones et al. 
(|2003l ). who defined a criterion of diffuse X-ray 
emission in excess of 0.5 x lO 42 ^^ 2 erg s _1 for 
such a classification. This arises from the process- 
ing of a group's IGM during a merger (or series 
of mergers), but the low total mass of most local 
compact groups suggests their potential well lacks 
the depth requ ired to heat the IGM to X-ra y de- 
tectable levels (|Mulchaev fc Zabludofilll998l ). Us- 
ing multiple mergers as a vehicle toward a fossil 
group end-state maps one path of galaxy evolution 
from the 'blue cloud' of star-forming disk galaxies 
to the 're d sequence' of qu iescent bulge-dominated 
galaxies (|Bell et al.ll2004h . 

Fossil group formation may provide an analogy 



2 



to cluster centers or sub-clumps where the buildup 
of cD galaxies occurs. If this turns out to be valid, 
the study of compact groups could also help illumi- 
nate morphological transformations in the inner- 
most cores of clusters. Exploring these different 
scenarios may prove fruitful for our understanding 
of galaxy evolution and the buildup of stellar mass 
in the universe. Making meaningful progress in 
this area requires detailed multi- wavelength stud- 
ies in order to map the range of physical processes 
affecting galaxies that are found in these environ- 
ments, determine their histories, and project their 
evolution. 

A consistent treatment of a large sample of 
HCGs is therefore i n order. In th i s wor k we 
continue the series of iGallagher et all (|2010h and 
Konstantopoulos et al. ( 2010l ) and provide a com- 
prehensive, multi-wavelength study of HCG 59. 
We will look at the current state of the group 
through its star formation and nuclear activity; 
investigate its past through the star cluster pop- 
ulations; try to unravel the history of mergers in 
the group; examine its dwarf galaxy system; and 
place it in the context of HCGs in general. 

The core of HCG 59 consists of four giant 
galaxies, a typical number for HCGs in general. 
The group lies at a distance of 60 Mpc, based 
on a recession veloc ity of vr — 4047 kms -1 
( Hickson et al.l Il992 . corrected to the reference 
frame defined by the 3K Microwave Background) 
and Hq = 73 kms _1 Mpc _1 . Three of the galax- 
ies, A (type Sa), B (E0), and C (Sc), have 
seemingly undisturbed morphologies, and the 
fourth (D, Im) is an unusually large irregular with 
a normal, peaked light profile. The total stel- 
lar mass of the group is M£ OT = 3.14 x 10 10 Mp, 



(from the 2MASS K s -ba,nd luminosities; lTzanavaris et al 
20101 ). while the Hi mass of -Mm = 3.09 x 10 9 M 



is comparable to the value expected for the mor 
phological types and ste llar masses of the member 



galaxi es, according to iVerdes-Montenegro et al 
(|200l[ ). This is therefore a somewhat gas-rich 



compact group, given that HCGs typically con- 
tain only about a third of the Hi expected. On 
the other hand, the J07 scheme classifies the Hi 
content of the galaxy group as a Type II, i. e. in- 
termediate in gas content, according to its ratio of 
gas-to-dynamical mass of log(Mni)/ log(Md yn ) = 
0.81 ± 0.05. These classifications are based on dif- 
ferent criteria and the disparity can thus be rec- 



onciled. Table [T] summarizes some of the general 
characteristics of the four galaxies, while Table ?? 
presents some derived and literature values of the 
mass content and nuclear identifications in the 
four galaxies. 

This paper is organized in the following way: 
Section [2] presents the optical, IR, and X-ray 
datasets used throughout this work. Section[3]pro- 
vides a full account of the young and old star clus- 
ter populations, which we use as our prime diag- 
nostics of current star formation and ancient inter- 
actions. In Section |4] we discuss the main findings 
of this work. Finally, in Section [5] we summarize 
the work presented and offer ties to previous and 
future work in this series. 

2. Observations 

2.1. HST optical imaging 

The analysis presented in this paper is based 
largely on HST- ACS /WFC multi-band data. Im- 
ages were taken in the F435W, F606W and 
F814 W bands in two pointings to cover all known 
giant group members. We will refer to these filters 
as £435, Ve>06) -^814 (and the set as BVI) to denote 
the closest matches in the Johnson photometric 
system. The notation does not, however, imply 
a conversion between the two systems. The ob- 
servations were executed on 2007 February 24, as 
part of GO program 10787 (PI: Jane Charlton). 
The exposure times were 1710, 1230 and 1065 sec- 
onds in the BVI bands respectively. Three equal 
sub-exposures were taken with each filter with 
a three-point dither pattern (sub-pixel dither- 
ing). Images were reduced 'on the fly' to pro- 
duce combined, geometrically corrected, cosmic- 
ray cleaned images. For the analysis of point 
sources, we used the standard HST pipeline prod- 
ucts with a nominal pixel scale of 0'.'05 per pixel. 
For analysis o f the extended sources , we ran 
MultiDrizzle (jFruchter fc Sosevl [2009h with the 
pixel scale set to 0'.'03 per pixel to improve the 
spatial resolution. The absolute image astrometry 
was checked with the world coordinate system of 
the Two Micron All- Sky Survey catalog (2MASS; 
Skrutskie et al. 20061 ) by identifying four unsatu- 
rated point sources in common; the average offset 
was ~ 0.01" in RA and Dec. The four main galaxy 
/814-band light p rofiles were fit with Sersic profile s 
using GALFIT (jSersid Il968t IPeng et al.l l2010al ): 
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Table 1: Basic information on HCG 59 main members 



Identifier 








Coordinates 3, 


Type 


m 


vr 


References' 3 












(J2000) 


H89 c 


(mag) 


( kms ) 




A: 


IC 0737 




11 


48 


27.55 +12:43:38.7 


Sa 


14.82 (B) 


4109 


[1], [2], [3] 


B: 


IC 0736 




11 


48 


20.08 +12:42:59.5 


E0 d 


15.60 (B) 


4004 


[4], [2], [5] 


C: 


KUG 1145H 


-129 


11 


48 


32.44 +12:42:19.5 


Sc 


15.90 (g) 


4394 


[4], [2] 


D: 


KUG 1145- 


-130 


11 


48 


30.64 +12:43:47.8 


Im 


16.00 (g) 


3635 


[4], [2] 



"Coordinates arc the centroids from fitting the HST TgM-band images of each galaxy with Scrsic profiles. See § 12.11 for more 
details. 

llli lEvans et al.N201CD ; [2l: lde Vaucouleurs et aTl Jl99lf ): [3l: lHickson et al.Nl992t) ; [41: lYork et all j200Cf l; [5l: lFalco et"aT1 dl999h . 
iHickson et~all ||1989T ) 

d The RC3 designation for 59B is' SO?', and so there is some uncertainty as to its classification. 




Fig. 1. — HST BVI color-composite imaging of HCG 59. The red colors of galaxies A and B imply little 
star formation, while C and D show some signatures of nebular emission in green. Also visible is the newly 
catalogued dwarf galaxy I (Section l4"7Tj) . 
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Table 2: Masses and star formation rates of HCG 59 galaxies 



ID 


M* a 


Af H2 b 


SFR C 


sSFR d 


Nucleus 




(xlO 9 


M ) 


(M Q yr- 1 ) 


(xl0~ n yr" 1 ) 




A 


17.40 


10.2 


4.99 ±0.67* 


28.66 


C/AGN 


B 


8.29 


< 8.8 


0.02 ±0.01 


0.19 


C 


C 


3.03 


< 7.2 


0.16 ±0.03 


5.15 


Hn 


D 


2.67 


< 6.6 


0.48 ±0.04 


18.15 


Hn 



a,c,d Stellar masses, star formation rates (SFR) and specific SFRs (sSFR) are drawn from Tzanavaris ct al. | teOlCT) . The published 
stellar masses were off by a factor of 7.4; these values have been corrected for this error. 
6 From the CO observations of Vcrdcs-Montencgro ct al. (1998) 

e From Martinez et alj (|2010T ) : 'C stands for 'composite', i.e. one that falls at the Hll/AGN overlap region, as defined in 
iKewlevTt" alTifeoOST ) ; see Sections 14.41 14.51 for discussion on these designations. 
■^This value is heavily affected by the AGN in galaxy A, as will be elaborated in Section 12.61 



the best-fitting centroid positions are given in Ta- 
ble Hi 

We used the images, presented in Figure [TJ to 
characterize the optical morphology of the galax- 
ies, and to detect and photometer star clusters 
and cluster complexes. All reported magnitudes 
are in the Vega magnitude system. In Section [3j 
we present the analysis of these two scales (clus- 
ters and complexes) of the star formation hierar- 
chy and also distinguish between young massive 
clusters (YMCs) and globular clusters (GCs). 

2.2. Optical point source photometry 

We follo w the same rationale applied in our pre- 

vious work ( Gallagher et al.l201dt Konstantopoulos et al 



and use star clusters to infer the star forma- 
tion activity and history in each of the HCG 59 
galaxies. At the adopted distance to HCG 59 of 
60 Mpc, we expect some contamination by super- 
giant stars, which can have absolute V"-band mag - 
nitudes as bright as —8.5 ( Efremov et alJ fl986h . 
At this distance, one ACS pixel measures ~ 13 pc, 
(cf. the average star clus ter radius of ~ 4 pc; e. g. 
Scheepmaker et al1l2007l ). meaning that clusters 
are at most marginally resolved and can be consid- 
ered point sources for the purposes of selection and 
photometry. We select cluste r s usin g the method 
described in iGallaeher et"aH J201Clh : in brief, we 
perform the initial selection on median-divided im- 
ages, require selection in all three bands, and filter 
the resulting catalog using point spread function 
(PSF) photometry. Our PSF filtering applied the 
following criteria from the output of the ALLSTAR 



routine in IRAP§ X values below 3.0; a sharp- 
ness in the range [—2.0,2.0]; and a photometric 
error less than 0.3 mag. Aperture corrections are 
first measure d between 3 and 10 p ixels and then 
added to the Sirianni et al.1 ( 20051) corrections to 
infinity. Finally, foreground (Galactic) extinction 
with E(B — V) = 0.037 is accounted for using the 
standard Galacti c extinction law (a correction of 
A v ~ 0.12 mag; ISchlegel et allll998l) . 

In order to fortify the selection against stars, 
we apply a conservative absolute magnitude cut 
at My < —9 mag, which produces the high- 
confidence sample. We do, however, define a larger 
sample by relaxing the magnitude cut and apply- 
ing stricter PSF-fitting criteria to detect globular 
clusters, which are expected to be fainter and 
point-like. In order to minimize contamination 
from marginally resolved sources such as com- 
pact b ackground galaxies, we follow R eikuba et al 



([20051 ) and apply hyperbolic filters (starting nar- 
row for bright sources and widening for fainter 
sources) with a maximum cut-off at the above 
mentioned criteria of magnitude error, \ an d 
sharpness. We will refer to this as the extended 
sample. 

The application of these criteria assigns 240 
bright star cluster candidates (SCCs) to the high- 
confidence sample and 948 to the extended sample. 
Specifically, the numbers of detected SSCs (ex- 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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tended sample numbers in parentheses) in galax- 
ies A through D are 7 (29), 77 (213), 13 (63) and 
65 (217), with a further 78 (426) objects coinci- 
dent with what would be the intra-group medium. 
We will provide a full analysis of these cluster pop- 
ulations in Section [3~T1 

In order to test the completeness of the final list 
of SCCs, we used ADDSTAR to add 3000 artificial 
stars to the image (over the entire field, includ- 
ing the galaxies) in the apparent magnitude range 
24-28, i. e. absolute magnitudes of (-9.89, -5.89). 
Because the final catalogue only contains sources 
detected in all three filters, we include this effect 
by calculating completeness fractions based only 
on artificial stars dete cted in all three bands (e. g. 
IPaRocha et all l2002h . The limiting magnitudes 
for the 90% and 50% recovery rates are (26.56, 
27.25), (26.51, 27.19) and (26.47, 27.15) in the 
B435 , V606 an d ^814 bands respectively (after pho- 
tometric corrections are applied) . For the distance 
modulus used of 33.89 mag, m — 26.5 mag corre- 
sponds to M ~ —7.4 mag. 

Our assessment of the state of star formation 
in HCG 59 is not limited to star clusters. Star 
cluster complexes represent a larger scale of star 
formation, as the optically blended concentrations 
of gas, stars, and dust that make up small star- 
forming regions, and likely include groups of clus- 
ters. In contrast to star clusters, these can be re- 
solved to even greater distances than studied here, 
as the fractal distribution of gas about a galaxy 
gives rise to such structu res at all scales (e. g . as 
demonstrated for M33 bv lBastian et al. 2007 ). 



2.3. Globular Cluster Candidate Selection 

Globular clusters are also selected from the 
HST images. Since the process is tuned to the 
color distributions found in HCG 59, we provide a 
full account below. As contamination from super- 
giants is not a problem for objects with GC-like 
colors, we adopt a fainter magnitude limit to se- 
lect old GC candidates than we used for SCCs. 
We have chosen a cutoff at Veoe = 26, which cor- 
responds to My ~ —7.7 at our adopted distance 
modulus for HCG 59, or slightly more luminous 
than the expected peak in the globular cluster lu- 
minosity function (GCLF) at My ~ —7.4 (e.g. 
Ashman fc ZepJll998t Irlarrid l200lh . The major- 
ity of GC candidates brighter than this limit lie 
above the 90% photometric completeness level in 



all 3 filters. Assuming a Gaussian GC luminosity 
function with a peak at My = —1 A ± 0.2 and a 
dispersion a = 1.2 ± 0.2, our faint-end cutoff then 
samples 39 ± 8% of the entire GCLF. 

We have selected GC candidates (GCCs) ac- 
cording to the color-space distribution of Milky 
Way GCs. We de -redden e d the colors of globular 
clusters from the lHarrisI (1996) catalog by their 
listed E(B — V) values, and then defined a paral- 
lelogram based on the intrinsic (B — V) — (V — I) 
color distribution of the MW GCs. This parallelo- 
gram was then converted to the A CS filters using 



the 's ynthetic' transformations in ISirianni et al 
(2005). The color selection region adopted here 
is 0.10 mag wider i n (V — I) than that used in th e 
analysis of HCG 7 ( Konstantopoulos et al. 2010h . 



but st il l doe s not exceed the boundaries of the 
lHarrisI (Il996h MW GCs. To qu antify, 95 of 97 
MW GCs from the lHarrisI (|l996l ) catalog (those 
with BVI information) lie within this box. All 
point sources in the HCG 59 fields with 1 a er- 
ror bars that overlap our color selection region are 
considered GC candidates, and are plotted in Fig- 
ure El 

Due to the close (projected) proximity of the 
galaxies in the group, it is likely that the halo 
GCs in each system will appear superposed. In an 
attempt to quantify the GCCs in each galaxy, we 
use the relationship between the galactic mass and 
the radial ex t ent o f the GC systems in galaxies of 



Rhode et al.l(|2007l ). To compute the expected size 



of each halo, we have adopted the mass-to-light 
conversions in that work, although we stress the 
general conclusions we reach are not dependent 
on the detailed size of any given halo. As the pre- 
dicted masses of all of the group gal axies are just 



below the lowest mass galaxies in the iRhode et al 



(|2007l) sample, we adopted a radial extent of 15 
kpc (or 56" at the assumed distance to HCG 59) 
for each of galaxies A, B and C. Galaxy D is a 
lower luminosity system, and as seen in Figure [5J 
GCCs in this object already lie within the pro- 
jected halo of GCCs in galaxy A. Discussion of 
the individual GC systems in each group galaxy 
follows in Section l3~3l 



2.4. Background and Foreground Contam- 
ination 

Contamination in our color-selected sample 
of GC candidates is expected from a variety of 
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Fig. 2. — Positions of detected globular cluster candidates, marked on HST 7gi4 imaging. We have also 
marked the expected GC haloes with blue boundaries. 'Outer B' refers to the location of a possible excess of 
GCs, discussed in the text. This is symmetric to the 'A-B bridge' (with respect to B; see Fig. 2]) and could 
be a tidally induced redistribution of the GC population (see Sections 13.31 and 14.41) . 
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sources, including foreground Milky Way halo 
stars, reddened young clusters, and unresolved 
background galaxies. With the small number of 
GC candidates present in some of the HCG 59 
galaxies (discussed below), contamination can be 
significant. 

Predictions from Milky Way star count mod- 
els (the Besancon model of lRobin et alj|20 03l) sug- 
gest that only 3-4 foreground Milky Way stars will 
appear in the magnitude and color ranges for ex- 
pected GCs in each of our ACS fields. Determining 
the contamination from younger, reddened clus- 
ters is more difficult, particularly in the central re- 
gions of the late- type galaxies C and D, where such 
objects could be present. Unresolved background 
galaxies are not likely contributing in any signifi- 
cant way to the numbers of objects in our fields; 
analyses of the bac kground objects (with GC-like 
colors) in HCG 7 (|Konstantopoulos et ail l2010h 
showed that the predicted foreground Milky Way 
star counts were similar to the observed putative 
background contamination, leaving little room for 
background galaxies to contribute significantly. 

We also consider the Pirzkal et al.l ( 20051 ) anal- 
ysis of stars in the Hubble Ultra-Deep Field 
(HUDF). Within the range of colors shown in 
Figure |8j they found the main contaminant of 
'void sky' to be M-stars, however, with a Veoe- 
Jgi4 of ~ 2.0, they are too red to be considered 
in our analysis. All Main Sequence stars detected 
in the HUDF are too bright to be mistaken for 
star clusters by our detection algorithm. In fact, 
the only class of stellar object that can be found 
in the color-space occupied b y our cluster can- 
didates is white dwarfs, which IPirzkal et al. find 
to have a density of 1.1 ± 0.3 X 10" 2 pc -3 . The 
maximum Galactic volume covered by our two 
pointings is a cube of ~ 8000 px on a side, or 
~ 0.24 pc~ 3 , assuming a scal e height of 400 pc (the 
maximum height quoted by IPirzkal et al. ) . Such 
a volume might be expected to host ~ 2.6 x 10~ 3 
white dwarfs. Therefore, we consider this poten- 
tial source of contamination to be negligible. 

The dominant source of contamination to our 
GCC sample will actually be stars from the Sagit- 
tarius dwarf galaxy. Our ACS fields (at I ~ 254°, 
b ~ +69°) are superposed on a rather dense part 
of the tidally induced (and bifurcated) lea ding 
arm (stre am 'A' from Beloku rov et al. 2006) of 
Sgr (e.g., iMaiewski et al.1 120031; iBelokurov et al 



20061: iNewberg et ail l2007t lYannv et al.l l2009h : 

the presence of such streams is not accounted 
for in traditional Milky Way star count mod- 
els. To investigate the impact of Sgr Stream 
stellar populations, i n Figure El we have overlaid 
12 Gyr isochrones of iMarigo et al. (2008) with a 
range of metallicities expected for t he Sgr lead- 
ing a r m, [M/H] 



Chou et al 



2007; lYannv et al 



-1 ± 0.5 (e.g. 

20091) . at distances between 



26 and 36 kpc onto color-magnitude diagrams 
of the point sources in both of our ACS fields 
(ass uming d ~ 31 kpc , with a spread of ±5 
kpc; Newberg et al.1 2007t iNiederste-Ostholt et al ' 



2010l ; ICorrenti et aLll2010l ). From this, we see that 
stars just below the Sgr Stream main sequence 
turnoff do have colors and magnitudes similar to 
that of the brighter (T^o6 < 25) GC candidates in 
our study, making some contamination likely. 

To estimate the total contamination in our GCC 
sample, we assume those GC candidates that lie 
far outside the GC system halos (as shown in Fig- 
ure [2]) are instead contaminating sources. The 
one exception to this is a region (called 'outer B' 
in Figure [5]) that lies outside the GC system of 
galaxy B, opposite to the direction of galaxy A. 
We will return to this feature below. There are a 
total of 18 objects in 8.4 arcmin 2 , or a background 
surface density £& a cfc = 2.3 ± 0.5 arcmin -2 . This 
is much higher than the predicted surface den- 
sity of MW halo stars from the Besangon model 
(Smw ~ 0.4 arcmin" 2 ), indicating that Sgr lead- 
ing arm stars are the dominant foreground source 
of contamination in our sample. 

Of course, for the above analysis we are mak- 
ing the assumption that these contaminating ob- 
jects are not bona-fide 'intra-group' GCs that lie 
far outside the main galaxies of the group. To test 
this, we compare the Veoe luminosity function for 
the background source sample with the luminos- 
ity function of the large GC candidate population 
surrounding galaxy B. We show this in the right 
panel of Figure [3] The extrahalo luminosity func- 
tion does not show a sharp rise with increasing 
magnitude as expected of a GC luminosity func- 
tion and exhibited by the GCCs in galaxy B. Al- 
though a definitive comparison is not possible with 
so few objects in the IGM area, this is consistent 
with these IGM objects being contaminants and 
not a part of a diffuse population of intragroup 
GCs. Thus we adopt the surface density above 



s 



HCG 59 




Fig. 3. — Diagnostic diagrams related to the GC popul ation. On the left we show color-magnitude diagrams 
of all detected GCCs. The solid and dashed lines show iMarigo et aL (2008) evolutionary tracks for 12 Gyr- 
old main sequence stars, with Z — —0.6 and —1.3 respectively. The lines are double to bracket the distances 
appropriate for the Sagittarius dwarf, with its leading spiral arm superposed along this line of sight. Their 
distribution in the CMD dispels our worries of significant contamination. The right panel shows the lumi- 
nosity function of globular clusters candidates (filled circles) in galaxy B and sources considered 'background' 
(open circles and dashed lines). The clear discrepancy supports the quality of our GC selection. The LF of 
'outer B' sources is designated with a solid line. 
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as that of the 'background' in the analyses that 
follow. 

2.5. Las Campanas wide-field imaging: 
low surface brightness light and dwarf 
galaxies 

We extend the coverage of the HST ob- 
servations through wide-field imaging with the 
Las Campanas Observatory (LCO) 2.5-meter tele- 
scope. We took B- and i?-band images of a 25' 
diameter around the group with the Wide Field 
Reimaging CCD Camera (WFCCD). The data 
were obtained on 2007 July 07 as part of an imag- 
ing campaign that covers all 12 HCGs in the J07 
sample. The B and R filter exposure times were 
300 s and 600 s, respectively. 

These images allow for the detection of low sur- 
face brightness features, such as the signatures of 
past interactions, over a large area. We present 
this analysis in Figure [4] where we stacked the B 
and R images and applied a Gaussian smoothing 
filter to the result. This image shows only fea- 
tures and i?-band contours that register at least 
3 cr above the background. We find two faint fea- 
tures, a 'bridge' that appears to connect galaxies 
A and B and an arc extending from B toward a 
compact structure to its north-west (we will later 
refer to this as the 'B-I arc'). There is another 
compact, extended source in the space between 
galaxies C and D. 

The original purpose of the LCO observing 
program was to prepare a sample of dwarf galaxy 
candidates for spectroscopic follow-up. Though 
our redshift survey has yet to cover HCG 59, it is 
covered by the S loan Digital Sky Survey (SDSS; 



York et al. 2000): a spectroscopic search sweep- 



ing a radius of 30 arcminutes around the nom- 
inal center of the group (the geometric center 
of the region enveloping the four known mem- 
bers) yields seven spectra with redshifts in the 
range 0.01-0.02: galaxies C and D and five com- 
pact galaxies. We therefore consider the mem- 
bership of SDSS galaxies J114817.89+124333.1 
and Jl 14813.50+123919.2, which are covered by 
our wide-field imaging and J114930.72+124037.5, 
J114940.11+122338.6 and J114912. 21+123753.8 
which lie at projected distances greater than 13 ar- 
cminutes from the group center. The first of these 
galaxies is also present in the HST imaging, but 
lies partly in the ACS chip gap. In Section 14+1 we 



attempt to determine whether these are HCG 59 
members through a phase-space analysis. 

The|Hickson| fjl982h naming convention assigns 
letters in order of brightness. Since our imaging 
does not cover all five dwarf candidates, we used 
the SDSS r-band photometry to consistently clas- 
sify the galaxies as HCG 59 F through J. We have 
omitted the letter E, as it was assigned in the 
original catalog to a background galaxy. We at- 
tempted to measure stellar masses for these galax- 



ies us ing 2MASS iiT-band images (jSkrutskie et al 
120061 ). however, they are below the detection limit 
of that survey. Table [3] summarizes all of the in- 
formation presented in this section: measured and 
SDSS photometry, radial velocities, galaxy mor- 
phologies and projected distances from the group 
barycenter. The latter two properties will be dis- 
cussed in Section |4~T1 



2.6. Spitzer observations: infrared spec- 
tral energy distributions 

The optical imaging was complemented by 
Spitzer imaging in the mid-infrared (IRAC 3.6— 
8 fim and MIPS 24 fim observations) presented 
in J07 and shown in Figure [5] In addition to 
the Ray leigh- Jeans tail of stellar photospheric 
emission, the IRAC bands probe the presence 
of hot dust and polycyclic aromatic hydrocar- 
bons (PAHs), while the 24 fim. observations trace 
cooler thermal dust emission. The dust and PAH 
emission are both stimulated by star formation 
activity. The harder spectra of active galactic nu- 
clei typically destroy PAH molecules while heat- 
ing dust to hotter temperatures than found in 
galaxies with star formation alone. At low AGN 
luminosities, the IR SEDs are often ambiguous 
(par ticularly in the prese nce of star formation; 
e.g- lGallagher et al.ll2008h . 

The Spitzer images were combined with JHK : 
observations from 2MASS (jSkrutskie et al.ll200 
to plot the IR spectral energy distribution (SED) 
of each galaxy (following J07), presented in the 
frequency-space p lot of Figure[6] We have used the 
Silva et al . (1998) templates for galaxies of vari- 
ous morphological types. These map the SED of 
different galaxies as the sum of starlight and gas 
and dust emission from star formation and inter- 
stellar cirrus. We calculate the spectral index of 
the SED within the IRAC bands t hrough a simple 
power-law fit. This was defined bv lOalla gher et al 
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Image (B + R) : 3 - 100 cr 
Contours (R) : 3 - 10 a 



Fig. 4. — Co-added, smoothed, background-suppressed B + R image of HCG 59 from the LCO duPont 
telescope. The image is aligned to the world coordinate system. The contours mark i?-band fluxes between 
3 and 10 a above the background (as marked by the circled numbers), while the smoothed B + R flux shows 
material brighter than 3 a above the background level. The asterisk label denotes the foreground star to the 
east of galaxy D, while the double asterisk marks the Hn regions west of B. The horizontal banding north 
of the main group is an artifact in the i?-band image. We detect a streak of material apparently connecting 
galaxies A and B in our field of view. We argue in Section l4~3l that this is a bona-fide tidal bridge, although 
its low surface brightness prohibits an in-depth photometric study. A faint arc of light to the west of galaxy B 
- containing the H II regions mentioned above (see § 13. ip - appears to connect it with the dwarf galaxy I. 
We discuss the importance of these features in Section [ 
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Table 3: Dwarf galaxies in HCG 59 



ID 


SDSS ID 


(Plate, MJD, Fiber) 


Morphology 






M r c 




dec 












(mag) 




( km s -1 ) 


(kpc) 


F 


J114930. 72+124037.5 


(1609, 53142, 464) 


Sd 


15.83 


0.09 


-18.15 


3984 


290 


G 


J114940. 11+122338.6 


(1609, 53142, 465) 


dir 


17.34 


0.08 


-16.63 


3195 


471 


H 


J114912. 21+123753.8 


(1608, 53138, 625) 


dir 


17.39 


0.08 


-16.58 


4011 


225 


I 


J114817.89+124333.1 


(1608, 53138, 586) 


dim 


17.43 


0.10 


-16.56 


4038 


32 


J 


J114813. 50+123919. 2 


(1608, 53138, 593) 


dE 


17.69 


0.10 


-16.30 


3942 


83 



"SDSS model r magnitudes. These values have not been corrected for Galactic extinction. 
^Galactic extinction and recession velocities are from SDSS DR7. 

c Absolute r magnitudes assuming a distance modulus of 33.89 and correcting for Galactic extinction. 





QJQQ ^^^^^ \ 



Fig. 5. — Spitzer 3.6-8 /j,m images of HCG 59. We present these individually to contrast the features evident 
in different bands and in a combined image to view the entire mid-IR picture. Galaxy A is extremely bright in 
the hot dust/PAH-dominated 5.8 and 8 fim bands. Galaxy B is faint in those bands, owing to little emission 
from heated dust. C is consistently bright in the four bands, given its ongoing star formation. D appears 
clumpy at long wavelengths, indicating patchy dust, heated by forming stars. In the 3.5 and 4.6/zm images, 
we also see faint, diffuse emission in the region between A and B, which we study in Sections 12.51 and 14.31 
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(|2008l ) as aiRAC, and it serves as a measure of star 
formation activity. In logarithmic frequency units, 
the flux difference from 8 to 4.5 fim leads to a posi- 
tive gradient in quiescent environments, while star 
formation registers as a negative slope. For HCGs, 
the steepness of the s lope is sensitive to the s pecific 



star formation rate ( Tzanavaris et al. 2010l ) 



In brief, the SEDs of galaxies B and C follow 
their morphological types of E/SO and Sc, while 
the irregular nature of D does not allow for a 
template to be assigned (we use Sc in Figure [6]). 
Galaxy A, nominally an Sa, shows strong excess 
light at wavelengths associated with PAH and/or 
hot dust emission (A > 5.8 /im). If the flux in 
this region is associated with star formation, an 
Sc template would be more appropriate. This is, 
however, inconsistent with the optical morphol- 
ogy of A. Furthermore, the shape of this emission 
is also consistent with the AGN dust b ump (e.g., 
Elvis et al.iri994lGallagher et al.ll2008l) . so in Sec- 



tion 14.51 we will examine the nuclear activity of 
this galaxy. 

The SED-fitting process could not be car- 
ried out for the new dwarf galaxies catalogued 
in Section 12.51 due to the limited spatial cover- 
age of the Spitzer imaging and the faintness of 
the dwarfs (they were not detected in 2MASS). 
As of July 2011, the Wide-field Infrared Survey 
Explorer (WISE) photometric catalogue did not 
cover HCG 59. 



2.7. The high energy picture: 
ACIS observations 



Chandra- 



The dataset is completed by Chandra-ACIS 
data in the 0.5-8.0 keV range. Results reported 
here are drawn from the work by Tzanavaris et al. (2011 
in prep.). HCG 59 was observed by Chandra be- 
tween 2008-04-12 and 2008-04-13 at the aim point 
of the back-illuminated S3 CCD of ACIS in very 
faint mode with an exposure time of 39 ks (ob- 
servation ID 9406, sequence number 800743, PI 
S. Gallagher). The data were processed using 
standard Chandra X-ray Center aspect solution 
and grade filtering, from which the level 2 events 
file was produced. Figure [7] shows an adaptively 
smoothed 3-band X-ray image, with optical (Isia) 
contours overplotted for comparison. 

Wavdetect, the CIA0 4 . 1 . wavelet detection 
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Fig. 6. — Near-to-mid IR SEDs for the primary 
HCG 59 members. The photometric data, shown 
as diamonds, are drawn from 2MASS (,JHK S ) and 
Spitzer MIPS/IRAC 3.5-24 ^m (from J07). We 
annotate each panel with the morphological type 
from the literature (top row), as well as the type 
chosen to most accurately represent the data (bot- 
tom row). The bracketed number quotes the age 
in Gyr of the template for the elliptical, or the 
inclination in degrees of the spir als. We also cite 



the v alue of the ckirac diagnostic ([Gallagher et al 



20081 ). which is fit on the wavelength range indi- 
cated by the dashed line. The solid curve s repre- 
sent GRASIL templates ( Silva et al. 19981) appro- 
priate for the morphological type of each galaxy 
(see text). Galaxies B and C are well represented 
by their nominal templates, while A and D de- 
viate: in the case of galaxy D this is due to its 
irregular nature, while we consider the nuclear ac- 
tivity of A to be the source of the excess mid-IR 
emission (see Section l4~5l for details). 



3 http:/ /cxc. harvard.edu/ciao 
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Fig. 7. — Adaptively smoothed three-band X-ray image of HCG 59, with overplotted /si4-band contours. 
Red, green and blue colors represent soft (0.5 — 2 keV), medium (1 — 4 keV) and hard (2-8 keV) X-ray 
emission. Galaxy C (Sc) has no detectable X-ray emission above the background level (the source at the 
limit of its optical isophotes is a likely background galaxy), while D (Im) shows some soft diffuse emission. 
The point source detected in the nuclear region of galaxy A (Sa) is detected in the full band with Lx (0.5-8.0) = 
1.1 x 10 40 erg s — 1 . This is surrounded by diffuse emission. We resolve two point sources in the central region 
of Galaxy B (E/SO) with Lx(o.5-8.o) = (1-4, 1.7) x 10 39 erg s _1 . Section 0~5] and Figure IT71 provide details. 
No diffuse emission is detected from the intra-group medium. 
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tool ( Freeman et al 1 l2002h was used in the soft 
(0.5-2.0 keV), hard (2.0-8.0 keV) and full (0.5- 
8.0 keV) bands to detect candidate point sources. 
The 1024x1024 S3 chip field was searched with 
wavdetect at the 10~ 5 false-probability thresh- 
old. Wavelet scales used were 1, 1.414, 2, 2.828, 
4, 5.657 and 8.0 pixels. Source lists produced by 
wavdetect for each band were matched against 
each other by means of custom-made scripts (K. D. 
Kuntz, priv. comm.) to calculate a unique posi- 
tion for each candidate point source, taking into 
account the varying size of the Chandra PSF 
across the S3 CCD. 

Point source photometry was carried out for 
the objects in the s ource list using ACIS Extract^ 
(|Broos et all l2010h which takes into account the 
varying Chandra PSF accross the CCD. Pois- 
son ±1(7 errors on net counts we re calcul a ted b y 
means of the approximations of ICehrelsl (|l98fih . 
Sources with measured net counts smaller than 
the 2 a error were flagged as non-detections. Note 
that this method produces very similar results to 
choosing a binomi al probability th reshold of 0.004 
in ACIS Extract (|Xue et al.ll201ll ). 

As sources have too few counts f or reliable spec- 
tral fit ting, we apply the method of lGallagher et al. 
(2005) to obtain a rough estimate of the spectral 
shape by using hardness ratios, defined as HR 
= (H - S)/(H + S), where H and S represent 
the counts in the hard and soft bands, respec- 
tively. Briefl y, we use the X-ray spectral modeling 
tool XSPEC (|Arnaudlll996h . version 12.5.0, to sim- 
ulate the instrumental response and transform ob- 
served HR values into an effective power law index 
T (where fx oc E~ T photons cm~ 2 s _1 keV -1 ), 
and also obtain associated X-ray fluxes and lu- 
minosities. This modeling includes neutral ab- 
sorption from the Ga lactic _Vh of 2.6 x 10 20 cm~ 2 
(|Kalberla et al.ll2005l ). 

We estimate flux limits of f x > 2.7 x 10~ 16 

erg cm- 2 s" 1 (0.5 - 2.0 keV) and fx > 1.6 x 

10~ 15 erg cm~ 2 s _1 (2.0 - 8.0 keV), corresponding 

l 

l 



to luminosity limits of Lx = 1.3 x 10 38 erg s" 



(0.5 - 2.0 keV) and L x = 7.7 x 10 38 erg s 
(2.0 - 8.0 keV). Assuming these limits, we u se the 



log N - log S relation of ICappelluti et al.l 12007) 



to estimate the number of background sources we 
would expect to detect in the HCG 59 field over 



the 1.96 x 10~ 2 square degree area of the S3 chip. 
This number is ~ 25 and ~ 20 in the soft and hard 
bands, respectively, with a ~ 20% uncertainty. In 
the much smaller area (3.66 x 10~ 4 square degrees) 
covered by our galaxies, we expect < 1 background 
source in each band. 

We detect a total of 40 sources in the soft band 
and 33 sources in the hard band over the ACIS S3 
field. We thus expect about 15 soft and 13 hard 
sources to be point sources associated with HCG 
59. We find that 9 soft and 3 hard sources are 
located inside the boundaries of the MIR-based 
HCG 59 galaxy regions of J07. We note that 11 
sources that have only hard-band emission are lo- 
cated far from the HCG 59 galaxies and are thus 
likely background AGN. 

In the central regions of two group galaxies, 
there are three notable X-ray point sources de- 
tected with high significance. As these point 
sources are all with 5' of the Chandra optical axis, 
the X-ray positions are the ACIS Extract "mean 
positions of events within the extraction regions" 
□ with intrinsic positional uncertainties of a few 
tenths of an arcsecond. There is some additional 
uncertainty from matching the absolute reference 
frames of Chandra and HST, but this is expected 
to be small as both are consistent with 2MASS at 
the - 0.1" level. 

The first source, in galaxy A, has a full-band 
luminosity of Ix(o.5-8.o) = 1-1 x 10 4 ° er S s 1 an d 
an estimated T = 1.3 ± 0.3. The X-ray posi- 
tion of this source is 0.7" from our quoted opti- 
cal position (Table 1). This isolated point source 
in the nuclear region of galaxy A has a luminos- 
ity that is consistent wit h known low-luminosity 
AGN (e.g jHo et afl200ll) and significantly higher 
than individual, luminous X-ray binaries. The 
two point sources found in the central region of 
galaxy B have X-ray positions 0.2" and 1.3" from 
the HST i-band galaxy centroid position (Ta- 
ble 1). Unfortunately, these sources have fewer 
than 10 counts in each band, precluding even a 
rough r estimate. Their full-band luminosities are 
£x(0.5-8.o) = (1-7, 1.4) xlO 39 ergs" 1 , respectively. 

The Chandra images are also sensitive to diffuse 
emission from MK degree gas. As can be seen in 
our adaptively smoothed image, some soft, diffuse 



*http: / /www. astro. psu.edu/xray/docs/TARA/AE.html 



5 See Section 5.3 of the ACIS Extract User Manual - 

http://www2.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
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emission is detected in galaxy D (Im) , likely asso- 
ciated with star formation, as well as in galaxy A. 
In both cases, the diffuse emission covers an area 
several times the size of the Chandra PSF at that 
location. We obtain an upper limit on the IGM 
surface brightness as follows. We calculate the 
count rate in a source-free region between group 
galaxies, and estimate the corresponding flux for 
kT = 0.5 keV thermal emission using PIMMS0 
We thus estimate the IGM surface brightness to 
be < 7.3 x 10 -17 erg cm~ 2 s _1 arcsec -2 . Finally, 
we note that dwarf galaxy I is not coincident with 
any detected X-ray sources; none of the other new 
dwarf members of HCG 59 are within the Chandra 
field of view. 

We will discuss the implications of these obser- 
vations further in Section 4. 

3. The young and old star cluster popu- 
lations: star formation over a Hubble 
time 

3.1. Young star clusters and the past ~ 
Gyr of star formation 

The population of star clusters is representa- 
tive of star formation as a whole in any system 
~ 2010h . 



(e. g. iBressert et al 



They are form ed 



masse after large events (e.g ., Trancho et al. 2007t 



Konstantopoulos et al.ll2008 , 20091 Konstantopoulod 
20091 : iBastian et all 20091 ) and at a slower pace 
at all times when a galaxy is forming stars 
(|Lada fe Ladall2003h . The extreme brightness of 
young clusters makes them detectable to large 
distances and therefore a reliable tracer of the 
star formation history of their host galaxy over a 
Gyr or so. Beyond that point in time they are 
referred to as intermediate age clusters and even- 
tually globular clusters. As a whole, the cluster 
population of a galaxy can reveal its star-forming 
history over a Hubble time. 

In this section, we analyze the cluster popu- 
lations of the four main galaxies in HCG 59 us- 
ing the HST PSF photometry described in Sec- 
tion 12.21 We use color-color (CC) and color- 
magnitude (CMD) diagrams to roughly age-date 
the clusters by comparing them to evolutionary 
tracks. The BVI filter combination of our HST 
images lacks coverage below ~ 4000 A, which is 



3 http: / /heasarc. nasa.gov/Tools/w3pimms. html 



crucial to breaking the age-reddening degeneracy 
(owing to the inclusion of the Balmer jump and 
near-UV continuum). However, it is still possible 
to infer the passing of intense bursts of star for- 
mation via the clumping of data-points along the 
evolutionary track, and unreddened young clusters 
are clearly evident. 

Figures [8] through [TT] show the high-confidence 
sample (i. e. clusters with My < — 9 mag; see Sec- 
tion 12. 2[) as solid dots, while the extended sam- 
ple is marked with open triangles. The solid and 
dashed red lines (r unning top to bottom) show 
Marigo et al. ( 20081 ) evolutionary tracks for simple 
stellar populations (SSP) of | Z Q and Z Q . This is 
sligh tly different from ou r previ ous work, where we 
used Bruzual fc Charlo tl (120031) models. We made 
this change because the lMarigo et al.l (|2008l ) mod- 
els seem to provide a good fit to both young clus- 
ters and globulars, unlike other model suites which 
focus on one part of the cluster population. The 
green lines that run more or less horizont ally show 
Starburst99 fSB99; lLeitherer et al.|[l999h tracks of 
the same metallicities. These also include nebular 
emission, which we expect to often be present dur- 
ing the first ~ 10 Myr of evolution. At this state 
the cluster is still surrounded by residual gas from 
the time of its formation ionized by UV photons. 
This short-lived phase ends when the first stars 
evolve and explode in supernovae that expel the 
gas. Clusters with colors redder than 0.8 in both 
axes are most likely GCs, although they might be 
highly reddened young clusters. 

When contrasting the extended and high- 
confidence samples, the former appears to spread 
more in color-space. This reflects a mass ef- 
fect intrinsic to star clusters, rather than in- 
dicating contamination. T wo recent studies, 
ISilva- Villa fc Larsenl ( 2011 ) and Popescu fc Hanson! 
(|2010h , independently reached the conclusion that 
lower mass (fainter) clusters, M < 10 4 M Q , of- 
ten exhibit deviations from the theoretical model 
tracks. In this mass regime, the underlying stellar 
initial mass functions (IMF) are under-sampled 
and stochastic effects dominate the overall light. 
Given that the IMF is populated randomly, it is 
physically equivalent to creating either one high 
mass star, or ~ 100 low mass stars (given the IMF 
slope). In contrast, a high-mass cluster will popu- 
late the IMF fully. Consequently, a population of 
high-mass clusters will have smaller intrinsic pho- 
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tometric dispersion. A cluster with lower mass 
will run out of material before the IMF becomes 
fully sampled, thus leading to the presence of gaps 
and spikes. Because of this effect, the presence of 
a high mass star in a low mass cluster will make 
the cluster appear to have a larger photometri- 
cally derived mass than one with only lower mass 
stars. 

We now treat the population of each galaxy 
individually. In the following paragraphs, all 
star and globular cluster candidates discussed are 
taken from the high-confidence sample. The dis- 
tinctions of young, intermediate, and globular are 
inferred from the locations of candidates along 
evolutionary tracks within the -B435-V6O6 vs. Veoe - 
7gi4 color-color space. 

Galaxy A (Sa) hosts a very small detectable 
population of five GCCs and two intermediate- 
age cluster candidates. The CMD shows them 
all to be consistent with having high masses, with 
log(M/M ) > 5. The lack of young clusters, com- 
bined with the high masses of the ones detected, 
indicates that galaxy A has stopped forming stars 
at a high rate. Because the mass-to-light ratio 
of a stellar population increases with age (i. e., 
the older a cluster becomes, the higher its mass 
must be for detection) , young clusters should dom- 
inate the cluster population when the SFR is high. 
The non-detection of any young clusters is fur- 
ther evidence that the UV+IR emission is domi- 
nated by AGN continuum, and the inferred SFR 
of 5 M Q yr _1 is severely overestimated. 

Galaxy B (E/SO) hosts mostly red clusters; 
the GCCs vastly outnumber young cluster candi- 
dates in this system. In fact, the four SCCs in the 
high-confidence sample do not strictly belong to 
galaxy B, but are found in the IGM: two nebular 
sources are located in a clump to the west of the 
galaxy, which we are treating as part of a stream 
that connects this galaxy to the dwarf HCG 591 
in our field of view. We will return to these H II 
regions in Section l4~3l The two non-nebular SCCs 
are distinct clusters in dwarf galaxy I. This means 
that the presence of patchy dust (as seen in the 
HST images of 59B) by itself is not indicative 
of star formation at a high enough level to pro- 
duce massive clusters, consistent with its SFR of 
0.02 M yr" 1 . 

Galaxy C (Sc) shows a continuous star forma- 
tion history, evident in the smooth distribution of 



datapoints along the evolutionary track. The pres- 
ence of nebular sources indicates some current star 
formation, while the low masses (« 10 4 M©) de- 
rived throughout the sample imply an overall low 
level of star formation over time. This is in ac- 
cord with the low value of 0.16 M Q yr _1 for the 
SFR of the galaxy. There is no pronounced GC 
population. 

Galaxy D (Im) is an unusually large irreg- 
ular galaxy. It shows a continuous star forma- 
tion history through to the present. The CMD 
shows a handful of SCCs with log(M/M ) > 5. 
At 0.48 M yr~l, this galaxy has the largest SFR 
of those with young clusters. There is no old com- 
ponent in the cluster distribution, no evident halo 
of GCCs, again possibly due to the low mass of the 
system which implies a small globular cluster sys- 
tem. Furthermore, the extended sample does not 
reveal a tight, correlated color distribution charac- 
teristic of GCs. The youngest clusters lie at a typ- 
ical color-space distance of ~ 0.3 mag away from 
the nebular model track along the reddening vec- 
tor, indicating the presence of dusty star-forming 
regions. 

3.2. Star cluster complexes 

One step above star clusters in the star for- 
mation hierarchy is cluster complexes, large ag- 
glomerates of young stars, arranged in a fractal 
distribution that follows the collapse of the pro- 
genitor gas. These structures can be used to un- 
derstand the global star formation activity in a 
galaxy. They are found to be more compact at 
higher redshifts than in the local universe (with 
two to fi ve times higher mass s urface density, as 
found by lElmegreen et ah 20091) . although in one 
local interacting compact group, HCG 31, we find 
complexe s to be similar to tho se at intermediate 
redshifts (|Gallagher et al.ll201dh . We argue in this 
series of papers that compact groups might process 
gas more efficiently when interacting than most 
other environment s apart from the in fall regions 



of galaxy clusters (jWalker et al.ll2010h . HCG 31 



fits well in that context. 

Star-forming complexes are extended amor- 
phous regions, with dust scattered across their 
surfaces. Their boundaries were identified by eye 
and measured by contours down to a limiting sur- 
face brightness about 10cr above the background. 
We single out 30 star-forming complexes in this 
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Fig. 8. — B435— V<506 vs. V(306— ^814 colors of all SCCs in HCG 59, p lotted along with i Z Q (solid red line, 
bottom) and Z© (dashed red) model tracks by iMarigo et al. ( 2008 ). spanning an age range of 6 Myr to 



13 Gyr (indicated by filled stars). The green lines (solid and dashed for | Z© and Z Q ) show a Starburst99 
track that accounts for nebular emission, as suited to the very youngest clusters that have not yet expelled 
their natal gas. Data-points that lie to the left of the dotted green line are considered 'nebular', i. e. younger 
than ~ 10 Myr. The high-confidence sample of SCCs (those with My < —9) is denoted by solid dots, 
while open triangles show fainter candidates comprising the extended sample (all cluster candidates to the 
detection limit). Much of the spread in the extended sample colors is due to stochasticity in populating the 
stellar initial mass function within clusters, as discussed in the text. We also indicate an extinction vector 
of length ^606 = 1 mag. We will present the distribution of each galaxy separately in the next figure. 
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Fig. 9. — S435— V606 vs- V(306— ^814 colors of SCCs, divided by galaxy, following Figured] Typical error 
bars are indicated as crosses in the lower left part of each plot. The irregular galaxy D is the only one 
that displays pronounced current star formation activity, as deduced by the number of nebular sources - 
although more nebular sources could be hidden behind dust and gas in the highly inclined spiral galaxy C. 
Galaxies C and D are consistent with continuous star formation histories over the course of their histories. 
The absence of GCs in D might indicate the galaxy is young, although it could be due to the fading of its 
low-mass cluster population below our detection limit (which is a function of age). Galaxies A and B only 
appear to host globular clusters and a handful of intermediate- age and young clusters, indicating quiet star 
formation histories over the past Gyr or so. 
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Fig. 10. — Vqqq vs. Vfsoe— Is ia color- magn i tude diagrams, following the datapoint plotting conventions of 
Figure |5J We have provided iMarigo et al. I (<2008h model tracks for four different cluster masses, times two 
metallicities, | Z Q and Z (top and bottom) (solid lines), with the 'nebular' tracks plotted as dotted lines. 
As in Figure|9l the age increases from top left to bottom right. In both samples, we tentatively deduce typical 
cluster masses - with the caveat that mass and extinction are degenerate - of < 10 M Q , with an additional 
locus of possible 10 6 M clusters. This distribution is consistent with local cluster populations, given the 
imposed cut-off for the high-confidence sample at My < —9 (and therefore the lower limit of ~ 10 4 M Q ). 
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Fig. 11. — V606 vs. Vqo6~ hu color-magnitude diagrams for each galaxy individually, following the plotting 
conventions of FigureQj)] Galaxy B is notable for the large population of apparently old, massive (~ 10 6 Mg) 
clusters, while galaxy D has shows a sizable population of 10 4 M Q clusters consistent with a predominantly 
young and intermediate-aged system. 
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compact group, three in galaxy C, and 27 in D. 
We find no complexes in the E/SO galaxy B or A, 
the evolved (Sa) spiral. The sizes of these com- 
plexes are com parable to those in lo c al sys tems 
(Figure QH cf. Elmegreen et al.|[l99l 1 19961) and 
consistent w ith HCG 7, a relatively inact ive com- 



pact group ( Konstantopoulos et al. 2010l ). Their 



colors are indicative of star formation and nebu- 
lar emission, as expected of star-forming regions 
(Figure [T21 right). In all, the star-forming com- 
plexes in HCG 59 are comparable to their counter- 
parts in local star-forming galaxies. They follow 
very closely a relation between brightness and 
size, much more so than we found in our previ- 
ous study of HCG 7. We show this correlation in 
Figure [12] (right) where the line is a simple linear 
fit to the data. The linear fit implies a similar 
surface brightness for complexes across the group, 
unlike in HCG 7, where we found significantly 
more scatter. 

3.3. The ancient globular cluster systems 
of HCG 59 

Globular clusters represent the earliest eras of 
star formation in a given galaxy. Their color dis- 
tributions carry the imprints of interactions and 
mergers and thus may probe the history of their 
hosts over a very long timescale. 

In this section we provide a full analysis of the 
number and color distributions of the GC popula- 
tions in HCG 59 and complement the star forma- 
tion and interaction histories we began to explore 
in Section 13.11 This analysis is built around the 
plots of Figure [13] The top row has a color-color 
plot (left) and histogram of the (B — I)q distribu- 
tions (right) in each galaxy. (The color-magnitude 
diagram is shown in the left panel Figure [3]) Ow- 
ing to their projected proximity, the GC systems 
of galaxies A and D overlap, and so we cannot pro- 
vide separate analyses. However, given the rela- 
tive masses of the two galaxies and their evolution- 
ary stages, it is safe to assume that the population 
is dominated by A. 

We find a fairly large population of globular 
clusters across this compact group, the vast major- 
ity of which are found in and around B. Recall that 
given the fading of clusters with age, GCs need 
extremely high masses, M > 10 5 Mq, to be de- 
tected to such large distances. We also find some 
clusters at large radial distances from the center of 



this galaxy, including some found along the stellar 
stream that appears to connect galaxies A and B 
(see Section 14. 3|) and its projection on the far side 
of galaxy B. Galaxies A and C host small pop- 
ulations. To quantify these populations, we de- 
rive the specific frequency, Sn, a measurement of 
the number of clusters per unit galaxy luminos- 
ity, for each galaxy. First, we correct the observed 
number of GC candidates by the background cor- 
rection noted above, and then calculate the total 
number of GCs expected around each galaxy by 
first adopting a photometric completeness fraction 
of / = 0.9 ± 0.1 for objects with Veoe < 26, and 
correcting for the expected fraction of GCs that lie 
below this magnitude limit. The Sn values are 0.3, 
7.7 and 0.1 for galaxies A, B and C respectively, 
assuming that all bright GCCs consistent with the 
halo of A are actually bound to A. The measured 
and derived numbers are collected in Table [4] 

The size of the population in galaxy C is consis- 
tent with its Sc morphological type, while B has a 
tremendously rich system, about twice the number 
of GCCs expected. Conversely, the GC population 
of galaxy A is much poorer than expected, com- 
pared to t he values of Sn ~ 1 typically seen in 
Sa spirals ( Chandar et al. 2004). Many of the de- 
tected GCCs appear to lie in a ring just outside of 
the bulge of the galaxy. Comparing to other GC 
populations in HCGs, our specific frequency for 
HCG 59B of S N = 7.7 ±3.0 is larger than (yet still 
consistent within the uncertainties) those observed 
so far in other large elliptical galaxies in HCGs: 
global values of Sn = 3.6 ± 1.8 for HCG 22A from 



Da Rocha et al 



HCG 90C from 



( 20021) Fl and Sn = 4. 4 ± 1.3 for 



Barkhouse et al 



pOOll ). 

Inspection of the region marked 'outer B' in 
Figure 3 shows a possible excess of GC candi- 
dates over the background level described above 
in Section 2.4, perhaps tracing an intra-group stel- 
lar population. This is not unusual, as compact 
groups by the ir very nature are likely to promote 
interactions. ( White et al. 20031) detected diffuse 
intragroup emission in HCG 90 accounting for up 
to h alf the light of the group. Two f u rther stud- 
ies 



.. . ,Da Rocha fc Mendes de Oliveiral (|2005h and 
Da Rocha et alJ ( 20081 ). found intra-group light in 



7 The lower Sjv ~ 1.7 ± 1.0 from Barkhouse et al. I d200ll) for 
HGC 22 A is due to their larger My for the galaxy — both 
studies find a very similar total number of globular clusters. 
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Fig. 12. — Color-color and size-luminosity (D — Mb) diagrams for star cluster complexes in HCG 59. The 
plotting of the evolutionary tracks follows the conventions of Figure [5] Complex symbols are coded to 
indicate where they are found: galaxy C (filled, yellow triangles) or D (filled, orange circles) The colors 
of the complexes are indicative of star formation, as expected. We find the complexes to follow a linear 
luminosity-size relation (see the upper left corner of the right-hand panel) very well. This implies a similar 
surface brightness for complexes across the group, unlike in HCG 7, where this relation is not obeyed. 




(B I). [M/H 



Fig. 13. — Left: Color-color diagram of all detected GCCs. T he dotted pa rallelogram on the left shows the 
selection box, based on the colors of GCs in the Milky Way ((Harris! 1996t h Center: {B — I)o distribution 
for all detected GCCs. The shaded area in the plot of the A/D distribution omits GCCs superposed over 
the body of D. The 'outer B' distribution covers all GCCs in pointing 2 that are at least 15 kpc from the 
centre of galaxy B. This region might suffer from some contamination, although it is consistent with observed 
distributions of GCs, which show a tendency for red clusters to be at large galactocentric distance as would 
be expected for a halo population. On the right, we relate the color distribution to that of metallicity and 
compare galaxies A and B directly. The dotted line in the GC metallicity distribution of galaxy B shows a 
best fit Gaussian, which we plot since a more complex distribution is not statistically justified. 
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another six groups: HCGs 15, 35, 51, 79, 88, and 
95. In HCG 59, a total of N = 15 ± 4 GCCs lie 
within the 2.7 arcmin 2 region, where we expect a 
background contribution of 6 ± 2 objects. Assum- 
ing only Poisson noise, this suggests a ~ 2a excess 
of objects in this part of HCG 59. This possible 
excess could be due to variations in the stellar den- 
sity in the Galactic halo or the Sgr Stream. How- 
ever, the luminosity function of this small number 
of sources in the 'outer B' region (shown as the 
solid line in the right panel of Fig. 3) is weighted 
towards the faint end, suggesting that at least 
some of these objects are indeed true GCCs lo- 
cated far (~ 25 — 50 kpc) from galaxy B. The un- 
usually large population of GCCs associated with 
B, the anomalously small population in A, and the 
possible population of GCs in the IGM (presum- 
ably stripped from a member galaxy) , all indicate 
that this compact group environment may have re- 
distributed GCs between member galaxies and/or 
to the IGM. 

In particular, possible interactions in the re- 
cent history of galaxy B might have redistributed 
its GC system. With this in mind, we compare 
the azimuthally averaged radial profile of the spa- 
tial distribution of the GC system with the surface 
brightness profile of the galax y, as derived through 
GALFIT (|Peng et al.ll2010bl ). The galaxy is best 
fit by a single-component Sersic profile of shape 
n = 3.1. Interestingly, we find the Sersic profile of 
the galaxy to provide a better description of the 
GC system than the best fitting power-law profile 
of index —1.3. This is contrary to the finding of 
power-law GC system distributions around loose 
group member NGC 6868 and HCG 22 A. If that 
is to be considered the norm for compact groups, 
then perhaps the recent interaction activity about 
galaxy B has changed the shape of its GC halo. 

The bottom panel of Figure [T3l shows the (B — 
I)o color distribution (converted to the John- 
son photometric system for direct comparison to 
Galactic globulars) of the clusters in each galaxy. 
The shaded area in the plot of A/D shows clusters 
that are clearly part of A, i. e. omits the ones that 
are projected upon the body of D. This does not 
alter the distribution, strengthening our assump- 
tion of a small population in D. Galaxy C also 
hosts a very small population, as expected due to 
its low mass. 

The color distribution can act as a proxy of 



metallicity for GCs and we take advantage of 
that to compare the two populous distributions 
of galaxies A and B. The color distribution of 
galaxy A seems fairly flat and the low numbers do 
not allow for a statistical treatment. Galaxy B, 
however, provides a large enough population to 
perform a test for bimoda l ity, u sing the KMM al- 
gorithm of lAshman et al.1 (|1994l ). This returns no 
evidence for a composite distribution in the metal- 
licity distribution. 

4. Discussion 

4.1. Extending the membership of HCG 59 



In Section [23] we introduced a search for dwarf 
galaxies in HCG 59, which we continue here. Since 
all objects we are considering here are covered 
by SDSS, we will not provide images and spec- 
tra here. More information can be obtained from 
the SDSS database using the plate IDs, Modified 
Julian Dates (MJDs), and fiber IDs given in Ta- 
ble E 

Regarding the morphologies of the new mem- 
bers, I, the candidate covered by our HST imag- 
ing, seems irregular, with a peaked light profile. 
This agrees with its spectrum, which shows clear 
emission lines and a continuum shape typical of a 
spiral. We classify it as dim. The rest of the galax- 
ies are not covered with high-resolution imaging, 
so we are more conservative in classifying them. 
F shows quite clear spiral structure in the SDSS 
images, and is therefore given an Sd type. We 
note that its star formation seems to be declin- 
ing, given that Ha is the only detectable emission 
line. G and H appear quite irregular and elon- 
gated, with spectra exhibiting blended emission 
and absorption. We assign them dlr types. Fi- 
nally, J shows a morphology closer to spherical 
and weak emission lines (although the S/N docs 
not allow for a confident determination); we as- 
sign it a dE type. 

Regarding their roles as group members, four 
of the five galaxies do not appear to be interact- 
ing with any other members, as might be expected 
from their locations far from the group core (Fig- 
ure HU left). It is only I that shows some evidence 
of an interaction with galaxy B, in the form of the 
'B-I arc' described in the Section [2~5l This stellar 
stream includes the large star-forming region we 
find to the west of galaxy B, which is part of the 
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Table 4: Observed and derived properties of the globular cluster systems in HCG 59 
Galaxy jV GCC a AWk b N v<2& c iV total d S N c 
A 18±4 6±1.0 13±6 33±21 0.3±0.2 

B 162±13 6±1.0 180±34 462±182 7.7±3.0 

C 7±3 6±1.0 1±3 3±8 0.1±0.3 

D 11±3 1±0.2 11±5 28±18 1.6±1.0 



"Number of detected GCCs 
''Background correction 

c Number of clusters brighter than V = 26 mag 
''Estimated total GC population 
e Specific frequency 



analysis of Section I3~T1 

In order to assess whether these galaxies belong 

to the group, we perform a pha se-space analysis, 

follow ing th e statistical studies of Mulchaev fc Zabludofj 
(|l998l ) and IZabludoff fc Mulchaev! (1200(1 ). This 
is shown in Figure [TJ1 right panel, where each 
datapoint represents a dwarf galaxy. The x-axis 
measures the distance of a galaxy from the group 
centroid, calculated as the mass-weighted average 
position of the four main members (A through 
D). The y-axis shows the offset of a galaxy's ra- 
dial velocity from the group mean, normalized to 
the core group (A through D) velocity dispersion 
of 314 kms- 1 . 

We find all galaxies to satisfy our member- 
ship constraints. Four of five spectroscopically 
detected galaxies lie within the boundary set by 
the four main members; galaxy G, which has the 
largest offsets in physical and velocity space, is 
moving with a radial velocity offset less than three 
times the group dispersion. We therefore con- 
sider all five galaxies under consideration here to 
be members of HCG 59, based on a strict 3 a- 
clip, as demonstrated in Figure [14] The inclu- 
sion of new members updates the velocity disper- 
sion of the group to 335 kms -1 . Based on this 
value, we derive a dynamical mass for HCG 59 
of M dyn = 2.8 x 10 13 M Q , a~ 10-fold increase 
with respect to the main members alone. This 
in fact changes the J07 evolutionary stage of the 
group from Type II (intermediate) to Type III 
(gas-poor), as it yields a ratio of H 1-to-dynamical 
mass of 0.71, with the caveat that the measured H I 
mass is likely underestimated because the dwarfs 
at large group radii are not included. 



It is interesting to find dwarf galaxies at large 
distances from the center of the group. This lack 
of barycentric clustering is also observed in the 
Local Group, where it is seen as a morphology- 
density relation: dwarf irregulars (dlr) ar e found 
at lar ger distances from the group center (Grebe! 
19991 ) than the quiescent dwarf spheroidals (dSph) 
and dwarf ellipticals (dE). This may indicate that 
some dlrs are galaxies experiencing their first in- 
falls to the group center. Such a situation could 
explain the relatively large velocity offset and the 
star- forming nature of galaxy G. 

The membership of galaxy G seems the most 
uncertain of the five galaxies discussed above, 
given the marginal agreement with the 3 a velocity 
cut, and the large projected barycentric distance. 
This is important, as its inclusion does affect the 
updated dynamical properties, due to the large 
change in group radius. We quantify this in Ta- 
ble El where we summarize the dynamical proper- 
ties of HCG 59. Those numbers show a change of 
mass by a factor of ~ 4 or ~ 10 with and without 
galaxy G, while the J07 type changes from II to III 
regardless of the inclusion of G. The velocity dis- 
persion is most affected: including G increases the 
value to 336 kms -1 from the original 314 km s" 1 , 
while excluding G significantly reduces the disper- 
sion to 208 kms" 1 , which is more consistent with 
the evolved state suggested by our analysis. 

4.2. The current state of star formation 

We have presented several diagnos tics of star 
forma tion activity across HCG 59. The lTzanavaris et al 
(|201dh SFRs of galaxies A through D are ~ 
[4.99,0.02,0.16,0.48] M Q yr" 1 . These are de- 
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Fig. 14. — Diagnostic diagrams of dwarf galaxy membership. The left panel shows the positions of the dwarfs 
and giant galaxies as blue dots and red stars respectively. We have labelled the dots with the appropriate 
lettering for HCG members. On the right we present a phase-space diagram that plots the offset of a 
galaxy's radial velocity from the established group mean (cr corc , the line of sight velocity dispersion of 
the giant galaxies) against its distance from the group barycenter (mass- weighted centroid). Four of five 
spectroscopically detected dwarfs lie within the boundaries set by the main galaxies and as such are confirmed 
as members. The fifth is moving at a relative radial velocity less than three times the group's velocity 
dispersion and is therefore also included as a member. The lower boundary of this 3 cr-clipping criterion is 
shown on the plot as the shaded area. 



Table 5: Updated dynamical properties of HCG 59 



Reference 


V 


a v 


Mdyn 


J07 type 




(kms 


- 1 ) 


(M ) 




Hickson fl982l 


4036 


314 


2.9 x 10 12 


II 


This work 


3924 


336 


2.8 x 10 13 


III 


This work, excluding G 


4015 


208 


1.1 x 10 13 


III 
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termined from the combination of UV and IR 
light assuming that all of the light emanates from 
star-forming regions. The presence of young star 
clusters in galaxies C and D shows that stars are 
forming at a fair pace, in accord with the star 
formation rates quoted in Table [TJ Overall, the 
SFRs are consistent with the infrared SEDs of 
these galaxies with the notable exception of A 
which is almost certainly strongly contaminated 
by AGN emission and shows no evidence of on- 
going star formation from the other evidence on 
hand. There are several soft X-ray point sources 
throughout the group, which are likely to probe 
compact stellar remnants local to HCG 59 galax- 
ies. There is also soft, diffuse X-ray emission con- 
fined to the galaxies, which probes ~ 10 6 K gas 
heated by star formation (stellar winds and SNe). 
The IR images do not show much that is surpris- 
ing: emission along the spiral arms of C and in 
the star-forming clumps embedded in D. Four of 
the five dwarf galaxies show star formation activ- 
ity, J being the exception. They are found to be 
star-forming based on either their bright emission 
lines or blue continua. Galaxies A and B, on the 
other hand, exhibit quiescent or even extinguished 
star formation. In the case of A, this is inferred 
by the absence of young star clusters. 

In all, the group does not appear to be un- 
dergoing a burst or any other event notable in 
terms of current star formation. With the ex- 
ception of the two most massive galaxies, the 
group is forming stars at a regular pace. This 
is also exemplified by the behavior of star clus- 
ter complexes across the group, which follow very 
closely after a brightness-radius relationship con- 



trast to HCGs 7 and 31 ( 


Konstantoooulos et al. 


2010t Gallaeher et alJl2010 


)• 



4.3. Signatures of interactions in the intra- 
group medium 

Major interaction or merger events are very 
often accompanied by bursts of star and cluster 
formation. The examination of star clusters in 
HCG 59 presented in Section I3TT1 did not show any 
evidence of such events in the last few Gyr. Given 
the high mass-detection limit for star clusters at 
this distance, they cannot be used to trace minor 
dynamical events. We therefore search for such 
evidence in the lowest surface brightness features 



detectable in our LCO images. In Section [231 we 
reported the detection of a low surface brightness 
stream of material in the projected area between 
galaxies A and B. This is visible at low surface 
brightness in our HST imaging and quite pro- 
nounced in the wide-field images from Las Cam- 
panas, at the > 3 a level. We also detected an arc 
of luminous material to the west of galaxy B, per- 
haps connecting it to compact galaxy I, which we 
discuss in Section |4~T1 

The low surface brightness and limited extent 
of the 'B-I arc' preclude precision photometry. We 
can therefore only pursue an in-depth analysis of 
the bridge between A and B. In order to derive 
the photometric properties of this feature, we first 
drew a color-map to look for an evident B — R 
gradient. Unfortunately, the stream is not bright 
enough to clearly dominate the image background. 
We thus conducted photometry of the area and 
the outskirts of the two galaxy that bracket it. 
We used large apertures of radius 75 pc (8 px) 
to reduce the background noise. The results are 
plotted in Figure 115) on the left panel we compare 
the measured photometry with the evolution of 
the B — R color, according to the iMarigo et al 



(|2008l ) model tracks of three metallicities. We find 
an intermediate value between the colors of the 
outskirts of the two galaxies that define this region 
(plotted as yellow dashed lines), cautioning that 
the emission in this region might be affected by 
the two galaxy light-envelopes to some extent. 

The origin of the bridge is not clear and its 
faintness makes it difficult to ascertain the domi- 
nant source of emission. It could consist either of 
stripped stars, or stars that formed in situ from 
stripped gas. If this is mixed stellar material from 
the two galaxies, we cannot study it in any more 
detail. We can, however, develop the in situ for- 
mation scenario further, by treating the bridge 
as a simple stellar population. In this case, the 
color-magnitude diagram of Figure [T5] (right), pro- 
vides an age estimate of about 1 Gyr, depending 
on metallicity. In addition, the CMD plotted in 
this figure provides an estimate of the stellar mass 
contained: with 10 7 M Q of ~ 1 Gyr old stars in 
each aperture, we extrapolate a mass in the order 
of - 10 s Mq, i. e. a density of - 100 M pc 2 . 
If the stars here are stripped from a galaxy, the 
overall mass will be higher, as the M/L of simple 
stellar populations increases with time. 
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Fig. 15. — Diagnostic plots of the tidal bridge between gala xies A and B . The age-color plot on the left 
shows the evolution of the B — R index over time (logr), for iMarigo et al. ( 20081 ) models of three different 
metallicities: I, \ and solar, shown as dashed green, dashed-dotted blue and solid red lines respectively. The 
range of colors in the bridge, as measured in five successive apertures and shown as a shaded region, lines 
up with the tracks at ages around 1 Gyr (depending on metallicity). The photometry of this faint feature 
might be affected by diffuse light from the galaxies, which brackets that of the bridge (denoted by dashed 
yellow lines, galaxy A is the bottom line). A secondary diagnostic is presented on the right, in the form of 
a color-magnitude diagram. The numbered dots denote log age in years, while the stars show the color and 
magnitude of each bridge aperture. If this is a simple stellar population (see text), then we can extrapolate 
a mass of ~ 10 8 M Q over the ~ 1 kpc 2 area of the bridge. This measurement can provide constraints for 
models of the evolution of tidally induced stellar structures. 
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4.4. HCG 59B as a merger remnant 

HCG 59B, the E/SO galaxy on the west side 
of the group, seems quite regular at first glance, 
however, a close inspection of the low-level light re- 
veals some interesting features. While ELLIPSE 
fitting shows an overall smooth isophotal struc- 
ture in B435 light, there is severe isophotal twist- 
ing in the central regions. This is spatially co- 
incident with several patches of extinction we de- 
tect in the HST images; they are most pronounced 
in -B435, observable in V606 an d hardly detectable 
in 1814) implying a thin column of dust. In the 
Spitzer bands, the fits indicate very symmetric 
structure in the 3.6 and 4.5 fjja bands, but the 5.8 
and 8.0 /im fits show a faint cusp some ~ 7 px, or 
~ 2.5 kpc from the r -1 / 4 surface brightness profile 
peak. 

The color-composite IRAC image (Figure [5]) 
shows hints of structure in galaxy B; however, 
an evolved elliptical/lenticular with near-zero star 
formation should present a smooth isophotal pro- 
file in all bands. Nonetheless, its SED (Figure |6]) 
shows a gradual decline, indicative of emission 
from stellar photospheres, rather than the heated 
dust associated with star formation. 

To investigate these irregularities further, we 
take advantage of the high resolution of the HST 
images. We construct pixel- by-pixel color maps 
of B, with the aim of tracing the exact location 
of the extinction patches. If the underlying stel- 
lar population is evolved to the same stage (i. e. 
an evolved SSP), then extinction will be the only 
source of discrepancies in color. There are, in fact, 
three possible sources of color variations in this fil- 
ter combination: (i) extinction across a similarly 
colored stellar population (mixed or coeval); (ii) 
spatially separated stellar populations of various 
ages and/or metallicities; or (iii) the presence of 
gas and star formation - i. e. H 11 region emission 
lines. 

The maps, shown in Figure [TBI cover the three 
possible filter combinations; for reference, the Veo6 
image is also shown. The B — I map is the 
one most sensitive to extinction. To quantify, 
the Cardelli et al. I (ll989l) extinction law assigns al- 
most twice the extinction in the -B435 as it does 
to the J 814 , A i35 /A 8 u = 1-85 (A 435 /A v = 1.13, 
cf. A S i4/Ay = 0.61). This is therefore the map 
we use to detect patches of low extinction, of 



order 0.3 mag, in three fingers extending ap- 
proximately eastwards from the north-south line 
through the nucleus. This faint structure is seen 
in all three colormaps, but not in the Veoe image. 
The F606W filter covers various emission lines, in- 
cludi ng Ha, [Nil], [S i ll , Hff , and [Om]. Interest- 
ingly, IMartmez etaD (|2010l ) found Ha, [Nil] and 
[SII] emission in the spectrum of B, at relative in- 
tensities consistent with a composite H 11 region 
plus AGN emission. From the concentrated, blue 
core of the V — I image, this line emission appears 
to be spatially coincident with the nucleus. 

Circumstantial evidence for a close encounter 
in the recent history of galaxy B is provided by 
the uneven distribution of its GC system. GCs 
normally form spherical haloes, however, here we 
find GCs at large radii, many concentrated along 
the stellar stream that seemingly connects galax- 
ies A and B, and its extension across the far side 
of galaxy B. Furthermore, the overabundance of 
GCs in galaxy B is matched by a severe dearth 
of clusters in A. Given the possibility that the 
two galaxies interacted ~ 1 Gyr ago, a scenario 
whereby GCs are transferred between the two sys- 
tems is not out of the question. It is unclear from 
a dynamical perspective why in the process of an 
interaction the GCs would flow from A, the more 
massive entity, to B. In a simple thought experi- 
ment, we move as many clusters from B to A as are 
required to level the Sn of A to the nominal value 
for an Sa. This still leaves an excess of GCCs in 
B relative to normal. However, the factor-of-two 
uncertainties involved in the determination of Sn 
do not rule out this scenario of GC 'swapping'. 

The X-ray map of this galaxy, as described in 
Section [2~71 reveals two distinct X-ray sources in 
the nuclear region. Unfortunately, with full-band 
luminosities of Lx = (1-4, 1.7) x 10 39 erg s , nei- 
ther is sufficiently luminous be identified as an 
unambiguous AGN - a possibility that the opti- 
cal spectroscopy of IMartmez et al" has indicated. 
Given the lack of ongoing star formation in the 
region, the sources are unlikely to be high mass 
X-ray binaries, though luminous low mass X-ray 
binaries (associated with older stellar populations) 
or groups of them unresolved at the distance of 
HCG 59 are plausible. In addition, due to the 
uncertainty in matching X-ray sources to optical 
imaging, we cannot confidently derive a one-to-one 
correlation between the optical clumps and these 
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Fig. 16. — HST color maps of the central region of HCG 59B. The maps are centered on the nucleus, defined 
as the peak of the 7si4 surface brightness profile, and aligned to the world coordinate system, with north 
to the top and east to the left. The bottom right map shows the Vqqq brightness. Contours show -B435, 
^oe, -^814) an d ^814 flux in the four panels (upper left, upper right, lower left, lower right). Each successive 
contour maps a difference in brightness of 0.5 mag. The B — I map on the bottom left is most sensitive to 
extinction and does indeed show patches of moderate reddening. The top row maps also show the extinction 
patches evident in the bottom row. In both the V — I and B — I color maps, the nucleus appears bluer than 
the surroundings; this is consistent with additional flux from line emission in the B435 (H/3) and Veoe (Ha) 
filters. The meaning of the observed colors is discussed in more detail in Section WM 
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sources, although the correlation is confirmed to 
within the X-ray positional uncertainties. We note 
however that, as discussed in Section 2.7, the nu- 
clear X-ray sources are spatially close to the opti- 
cal center of the galaxy. 

The excess of globular clusters and uneven dust 
distribution in the nuclear regions of B hint at 
some interaction in the more distant (> Gyr) 
past, but the lack of additional evidence for struc- 
tural disturbances limits our ability to infer more. 
We do detect a few young clusters in galaxy B 
and a non-zero (though low) SFR, and therefore 
some reservoir of cold gas is present. Accretion 
of a satellite galaxy is therefore a possibility. 
Furthermore, the unimodal GC color distribu- 
tion does not favor a gas-rich, major merger in 
th e past. This conclusion f ollows the paradigm 



of iMuratov fc Gnedinl (|2010l ) , who attribute the 
known bimodality of GC colors to late-epoch 
mergers. It is also a reasonable assumption that 
the implied interaction did not feature a major 
merger with a gas-rich system, as that would have 
enhanced the young and intermediate-age cluster 
populations. 

4.5. Nuclear activity in HCG 59A 



In Section [2761 we reported that the IR emission 
in galaxy A is consistent with being dominated by 
an AGN, rather than star formation. This is based 
on the disparity between the galaxy's morphologi- 
cal type of Sa and the high UV+IR SFR suggested 
by interpreting the emission as related to SF. Fur- 
thermore, the lack of young massive clusters is in- 
consistent with a SFR of 5 M Q yr" 1 . 

This hypothesis is supported further by the 
finding of a hard X-ray source in the nuclear re- 
gion of galaxy A, with Lx = 1.1 X 10 40 erg s" 1 as 
reported in Section 12.71 The spectros copic AGN 
survey of HCGs (|Martmez et al.ll2010h places the 
galaxy at the interfac e of the H II and AGN zones 
in the 'BPT' diagram (jBaldwin et al.lll98ll ). based 
on optical emission-line ratios. The 2"-wide slit 
they used encompasses a wide region (effective 
aperture of 0.58 kpc), therefore the signal is most 
likely diluted by circumnuclear and disk light. 
Visual inspection of the central region reveals 
asymmetric structure, resembling a second nuclear 
source of comparable /si4 luminosity to the nu- 
cleus. As in the previous section, we employ HST 
color maps to take advantage of the spatial reso- 



lution of ~ 12 pc per pixel. 

The maps, shown in Figure I17[ reveal a cone 
of blue light at (x 7 y) ~ (—0.1,0.1), with colors 
of B 435 - V eoe ~ 0.9, - hu — 0.6, and 

-B435 — -^814 — 1-6- The complexity of the central 
region inhibits easy interpretation, but one clue 
are the emission lines covered in the three bands. 
Vf306 covers H/3, [Om], Ha and [Nil], which can be 
associated with star formation and/or AGN activ- 
ity. In this scenario, the blue cone could stand out 
as a result of geometry, perhaps being located in 
a break in the dust distribution. A simple expla- 
nation could relate this feature to an unreddened 
line of sight through the inner spiral structure of 
A. In that case, however, we would expect to see 
bright, young star clusters, as they are known to 
shine through thick columns of dust, let alone rel- 
atively dust-free regions (e. g. Region B in M 82 ; 
Smith et al.ll2007t iKonstantopoulos et al.ll2008l) . 



A different interpretation can relate this struc- 
ture to an AGN. The lack of symmetry could 
suggest a small narrow-line region, photoion- 
ized by the AGN continuum, with projection ef- 
fects and obscuration hiding the cone on the far 
side. This would produce strong [O ill] and [N 11] 
emission, the p r esenc e of which was reported by 



Martinez et al 



<l2010h . This geometry is consis- 
tent with the inclination of galaxy A of no more 
than 30° (assuming the AGN and galaxy share 
the same inclination angle). 

Combining the pieces of evidence collected from 
the X-ray, optical and MIR emission, we propose 
that the nuclear emission in HCG 59A is domi- 
nated by a low-luminosity AGN with a photoion- 
ized narrow-line region. The onset of activity may 
be related to a possible encounter with galaxy B 
about 1 Gyr ago, as tentatively dated from the 
colors of the bridge connecting the two galaxies 
(Section O. 

5. Summary 

We have presented an analysis of HCG 59, 
a compact group comprising four main galax- 
ies and at least five newly discovered dwarfs at 
the M r < —15.0 mag level. Our results are 
based on multi-wavelength observations and con- 
tinue a series of papers that have followed two 
different approaches: on the one hand we have 
treated the overall properties of Hickson compact 
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Fig. 17. — HST color maps of the nuclear region of HCG 59A, following the plotting technique presented in 
Figure [16l North is to the top and east is to the left. The central peak appears elongated in the Veoe flux 
map (lower right) and shows a composite structure with two additional peaks nearby and a long tail. The 
tail could be an artifact introduced by the complex and heavy extinction in this region, perhaps caused by 
the inner spiral structure of HCG 59A. It could alternatively be associated with the nuclear activity we note 
in the X-ray image of this galaxy (Section 12. 7|) , as well as the significant excess in 8 and 24 /im emission, 
discussed in Section l2~6l A blue cone-shaped structure is most evident in the Veoe-Igu color map (top right). 
The nature of this blue cone is discussed in Section |4~51 
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groups rtJohnson et al.)l2007tlGallagher et al.ll2008t 



Tzanavaris et al.ll20ialWalker et al.ll2010t ): on the 
other hand, we have surrounded our HST obser- 
vations with a multi- wavelength dataset to pursue 
in-dept h investigations of i ndividual CGs, one at 
a time Jpalma et al.ll2002t iGallagher et alj|20ld: 
Konstantopoulos et al]l201of h 

Compared to HCGs 7 and 31, two compact 
groups previously studied in this series, HCG 59 
presents something of an intermediate step: where 
HCG 7 was found to be interacting solely in the 
dynamical sense (i. e. currently in the absence of 
direct hit encounters), HCG 59 shows evidence for 
stronger interactions in the recent past. There 
is evidence for star formation in the intragroup 
medium in the H n regions to the Northwest of 
galaxy B, in sharp contrast to HCG 7 where the 
star-formation associated with each galaxy was 
self-contained. In the cont ext of the evolutionary 



sequen ce we proposed in iKonstantopoulos et al 
(l2010h . it occupies a stage further along than 
HCG 7. It has begun building an IGM, as tes- 
tified by some amount of intra- group light, but 
has yet to build up a large elliptical fraction. Its 
classification as a relatively unevolved group (J07) 
is in accord with standard diagnostics such as its 
velocity dispersion of ~ 335 kms -1 and lack of 
diffuse, extended X-ray emission. 

Through the use of SDSS data, we associated 
five dwarf galaxies to HCG 59 for the first time. 
Their inclusion updated the velocity dispersion 
and dynamical mass of the group and changed 
its J07 evolutionary type to an evolved group of 
Type III (originally Type II, or intermediate). The 
star-forming nature of these dwarfs, their radial 
velocities, and distances from the group core seem 
to suggest that some may be infalling for the first 
time. The noted relation between morphological 
type and barycentric distance follows the one ob- 
served in the Local Group and may be considered 
- conversely - as an indication of the 'young' dy- 
namical state of HCG 59. The above information 
highlights the importance of studying the dwarf 
galaxy contingent of compact groups. 

Star formation is proceeding at a regular pace 
in this CG, certainly at a rate consistent with the 
morphologies of the member galaxies. The reg- 
ularity of star cluster complexes agrees with this 
image. The star cluster population does not show 
evidence of major, gas-rich, interactions in the 



past ~ Gyr and the IR SEDs are generally as ex- 
pected. One exception to this rule is galaxy D. 
Given its large size, the morphological regularity 
of its neighbors, and the lack of evidence of recent 
interactions, it is not clear why it has such an ir- 
regular structure and why it is forming stars at 
the rate that it is. 

Where the information is unclear for galaxy D, 
B is evidently in the midst of at least one dy- 
namical process. This probably started more than 
a Gyr ago given the lack of tidal features such 
as shells and tails co mmonly observed in such 
events in the optica l ( Schweizer fc Seitzer 19981: 



Iffibbard et all 120011 : iMullan et al.1 120111 1. It is 
likely physically associated with an arc of star 
formation that seemingly connects it to dwarf 
galaxy I (the arc that hosts the discovered ex- 
tragalactic Hn regions), and 59B lies at one end 
of a low surface brightness bridge of stellar ma- 
terial, which might physically connect it with A. 
We found this feature to be no older than ~ 1 Gyr 
by age-dating its stellar population. Perhaps most 
intriguingly, the globular cluster population of B 
is anomalously large, with significantly more glob- 
ular cluster candidates than A, despite its lower 
stellar mass. The origin of this discrepancy is un- 
clear, but hints at an additional event in the more 
distant past. 

In one sense, this group is typical of early or 
intermediate-stage HCGs. There are plenty of dy- 
namical processes at play, however they are all 
proceeding at a low level and centered around one 
or two objects. Galaxy B is at the focus of all 
such processes that our diagnostics can reach, like 
galaxy B in HCG 7 (also an early- type galaxy). 
In addition, both these groups may feature an ac- 
tive nucleus host galaxy (the dominant galaxy A 
in both groups, albeit the detection is not certain 
in HCG 7). 

We find the presence of a low-luminosity AGN 
in A to be likely given the spatial coincidence of 
the galaxy centroid with a ~ 10 40 erg s _1 X-ray 
point source. The inferred SFR from the IR+UV 
luminosity is not supported given the lack of young 
star clusters, and therefore likely is overestimated 
because of significant AGN contamination. If the 
A-B bridge constitutes a physical connection be- 
tween the two galaxies, a causal connection be- 
tween the interaction and the AGN is possible. In 
this case, we can constrain the timescale since the 
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onset to no more than one Gyr. 

In the introduction, we discussed HCGs as a 
potentially special environment in terms of galaxy 
evolution. The apparent duality of 'modes' in 
which HCG galaxies are found - either star- 
forming or quiescent - and the evident lack of an 
intermediate stage population are in accord with 
the mid-IR color-space 'gap' discussed in our pre- 



vious work (Johnso n et al . 2007; iTzanavaris et al 



20ld; IWalker et al.l l2010h . The impact of the 
group environment will be the topic of the 
next paper in this series. There we will treat 
HCGs 16, 22 and 42 with a goal of understand- 
ing the evolutionary processes at play in compact 
groups and relating their galaxies to those found 
at other levels of galaxy clustering. 
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